Polychromatic nanofocusing of surface plasmon polaritons 
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We introduce the concept of polychromatic plasmonics and suggest an broadband plasmonic lens 
for nanofocusing of surface plasmon polaritons. The lens employs a parabolically modulated metal- 
dielectric-metal structure. This plasmonic lens has a bandwidth of more than an optical octave thus 
opening new opportunities for broadband plasmonic applications. 
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The field of plasmonics experiences an explosive 
growth recently due to the ability of plasmonic compo- 
nents to confine light down to the nanoscale. Various 
plasmonic structures with miniaturization scales compa- 
rable to those of modern semiconductor electronics have 
been designed to realize different types of light waveguid- 
ing and control [H-Q- Despite the boom in research one 
can identify two major challenges for further applications 
of plasmonic devices: high propagation loss and strong 
wavelength dispersion. Different approaches to combat 
optical losses in metals have been explored, including in- 
corporation of gain [i| or nanofocusing 0, [Toj| . How- 
ever, the concept of broadband control and nanofocusing 
of polychromatic surface plasmon polaritons (SPPs) re- 
mains practically unattainable. 

In this Letter, we introduce the concept of polychro- 
matic plasmonics and demonstrate the functionalities of 
a broadband plasmonic lens based on a metal-dielectric- 
metal (MDM) structure (Fig. [I]). We utilize quadratic 
modulation of the thickness of the dielectric layer in 
transverse direction [Fig. GJb)] to produce a parabolic 
optical potential which is practically wavelength inde- 
pendent. We develop analytical descriptions and em- 
ploy numerical simulations to show its capability of 
three-dimensional subwavelength manipulations, includ- 
ing nanofocusing, self-collimation, and optical pendulum 
effect. The nanofocusing of our lens is demonstrated over 
a bandwidth exceeding an optical octave (> 500 nm) thus 
allowing for polychromatic plasmon focusing. 

The concept of polychromatic light propagation is well 
developed for dielectric structures, of which the polychro- 
matic dynamic localization is a noticeable example [Tl| . 
In these structures, curved waveguides are employed to 
produce an effective wavelength invariant optical poten- 
tial to compensate for the dispersion. 

In plasmonics, the engineering of the optical potentials 
for SPPs has been a key concept for designing plasmonic 
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lenses and other elements Jpl, |6|, [l2j. Two kinds of struc- 
tures with transverse 0, 0] or longitudinally varied ef- 
fective index [q|, [l3| are mainly used. In the structures 
with transverse index modulation the plasmonic poten- 
tials are usually heavily wavelength dependent [4-7] . On 
the other hand, in longitudinally modulated structures 
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FIG. 1: (Color online) (a) Flat and (b) parabolically modu- 
lated MDM. (c) Real part of the symmetric mode index of the 
flat structure vs. 1/h for various values of free space wave- 
lengths, (d) Real part of the effective permittivity (optical 
potential) of the structure (b), obtained using Eq. ([2]) (black) 
and numerically (red). 



light of different wavelengths can be focused asymptot- 
ically [l3| in the same taper, however light is fully lo- 
calized there and could not propagate beyond the focus 
point [sl, [l3( | . Recently a variation of a plasmonic lens 
was shown to focus broadband light based on diffraction 
[rH - However, the focusing was obtained for light in free 
space rather than for propagating SPPs, thus with low 
focusing resolution. The focusing of polychromatic SPPs 
still remains an unsolved challenge. 

To address this challenge, first we consider a flat MDM 
waveguide shown in Fig.QJa). We only consider the sym- 
metric modes of the structure (with respect to the mag- 
netic field distribution) because it is primarily excited by 
simple end-fire coupling, while the antisymmetric mode 
experiences a cut-off for the parameters of our work Q. 
The dielectric is chosen to be silica glass with permittiv- 
ity e = 2.25, and the metal is silver for which we use the 
Drude model, e m = 1 — uo^/(uo 2 + iujuj c ), where co>, uj Pi 
uj c are the angular frequency, plasma frequency and col- 
lision frequency, respectively. This model is proven to be 
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a good approximation for noble metals, including silver 
in the spectral range above 500 nm. For silver we take 
uj p = 1.37 x 10 16 rad/s and uj c /uj p = 0.002. Fig. HJc) 
shows the real part of the mode index, n e ff, versus the 
inverse thickness of the dielectric for three different wave- 
lengths. The dependence is linear, consistent with the 
theoretical approximation [Til, [l6| , 



n e ff = a/h + 6, 



(i) 



where both a and b are complex parameters and could be 
extracted from data fitting. Most importantly, the slope 
a of these curves is practically wavelength independent 
which allows for the design of broadband optical potential 
and polychromatic plasmon propagation. 

Using the effective refractive index in Eq. (pQ), it is 



possible to construct a parabolic optical potential [17] 
for SPPs. This is achieved by a MDM waveguide with 
one flat and one parabolically curved surfaces as shown 
in Fig.QJb). The thickness of the dielectric is h(x) = ho + 
x 2 /2Ro, with Ro as an effective radius (Ro ^> ho, \x\). In 
this waveguide, we obtain a parabolic optical potential 
under the condition of x 2 <C 2hoRo- 



£es = n e ff( x ) = n 2 (l-Q 2 x 2 ) 



(2) 



where no = a /ho + b, and Q = ^ 'a / '(noRoho) is the 
focusing strength. Fig. QJd) shows the results calculated 
by Eq. (j2j) and numerically for ho = 50 nm, A = 0.7 um, 
and Ro = 100 /im. 

Using the effective index method for the structure 
shown in Fig. QJb), we can express the vertical elec- 
tric field as E y (x,y,z) = A(x)B(x,y) exp(if3z) [l6l fl8|. 
When \x\ <C Ro, A{x) and B(x,y) could be decou- 
pled [16|. The expression for B(x,y) could be found in 
Refs. [J 16], while the equation for A(x) is 17, 19| 



d 2 A(x) 
dx 2 



+ [C 2 {x) - f3 2 ] A(x) = 0, 



(3) 



where Q(x) = kn e ^(x) and ("o = kno, k = 2tt/X. This is a 
harmonic oscillator equation with eigenmodes A m (x) — 
(^o2 m m!)-^ m (x/7 ?0 )exp(-x 2 /2^) [l7|,[l9|, where 
Vo = (Co^)~2 is the characteristic width of the plas- 
monic waveguide, H m is the Her mite polynomial and the 
effective indices of different modes, under the paraxial 
approximation, are n e ff (m) = [(o — \(^ m + 

We also calculate the eigenmodes of this struc- 
ture using commercial Mode Solutions (MS) software 
(Lumerical) with the parameters x G [—3, 3] /im, y G 
[—1, 1.05] /im, Ro = 100 /im, ho = 50 nm and A = 0.7 /im. 
The results of our analytical theory [Eq. (j3])] and the MS 
calculations are summarized in Fig. O Fig. [2fa,b) shows 
the real and imaginary parts of n e ff for different modes 
with the theoretical field distribution of \A m (x)B(x,y)\ 
[Fig. Efc-e)] and numerically calculated total field \E y \ 
[Fig. [2ff-h)], which show a good agreement for low or- 
der modes. For higher order modes, larger discrepancy 
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FIG. 2: (Color online) (a) Real and (b) imaginary parts of 
n e ff for different order modes. Theoretical (black circles) and 
numerical (red circles) results are shown, (c-h) Transverse 
mode profiles of the fundamental, second and tenth order 
modes, (c-e) Theoretically and (f-h) numerically calculated 
mode profiles (|i^|). 

appears in terms of both effective index and field distri- 
bution. This is consistent with the results in Fig.QJd), as 
the higher order modes spread out to the larger x values, 
where the assumption x 2 <C 2hoRo does not strictly hold 
and the potential is not exactly parabolic. 

Next we study the SPPs propagation in the structure. 
An incident beam could excite modes of different orders, 
which will then interfere with one another, producing dif- 
ferent intensity patterns inside the structure. As a result 
of this interference, a range of SPPs beam manipulations 
is possible, including focusing, self-collimation, and opti- 
cal pendulum effect [17|, [l9| . 

When A(x) can be decoupled from B(x, y) and B(x, y) 
changes slowly with x, the input beam expansions along 
x and y inside the structure are independent. Consid- 
ering that the beam distribution along y can be fully 
characterized by B(x,y), we study only the beam dy- 
namics in the x — z plane. The beam could be expanded 
into a complete set of orthogonal modes [l9[ C(x,z) = 
Em=o a m^mW exp(i/3 TO z), where a m is the expansion 
coefficient of the m-order mode. If the initial beam has 
a Gaussian distribution C(x,0) = tt~^w 2 exp[— (x — 
xo) 2 /2wq], where wo is the beam width, the beam inside 
the structure, under paraxial approximation, is [TR [l9|: 

C(x,z) = n~iq(z)~i exp j-^r [x - x cos(^z)] 2 | ^ 
x exp [i(k z + ^k p~ 1 x 2 )] P(x ), 

where rf = \q(z)\ 2 , p~ x = r]~ 1 dr]/dz, P(xq) = 
ex.p{i(r]o/wo) 2 [xxos'm(Qz) — -Xq s'm(2Qz)]} and q(z) = 
wocos(Q z ) + i{r]o/wo) sin(f^). In the lossless case with 
uo c = 0, the beam intensity can be expressed as: 
| C (x, z) | = 7r~irj~i exp j — [x — xq cos(Qz)] 2 / 2rj 2 j with 
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Wq = 340 nm w = 800 nm w = 480 nm w = 480 nm 
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FIG. 3: (Color online) (a-d) Theoretical field distribution 
without loss and (e-h) FDTD simulations with loss (shown 
\E\) in the x — z plane in the middle of the gap (y = ho/2 = 
25 nm). Plots (a, b, e, f) show typical effects of self-imaging 
with nanofocusing; (c, g) self-collimation; and (d, h) opti- 
cal pendulum effect. The corresponding characteristic beam 
width is 770 ~ 480 nm. 
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FIG. 4: (Color online) (a) Focusing position z — F\ for wave- 
lengths in the range 0.5 — 1.0 /mi at ho values of 20 nm, 30 nm 
and 50 nm. Solid lines: theoretical results; circles: FDTD 
simulations, (b-d) SPPs propagation for the three wave- 
lengths marked in (a). 



dynamic beam width varying along propagation: 

r] = yj Wq cos 2 (Qz) + Vo/ w o sin 2 (ftz) . (5) 

The plasmon polariton beam is thus oscillating periodi- 
cally along z. It is easy to prove that the maximum and 
minimum dynamic beam widths are constrained by [17| 
^max^min = Voi indicating that the beam is trapped in 
the parabolic optical potential. 

To confirm our theoretical analysis, we preform finite- 
difference-time domain (FDTD) simulations (Lumerical). 
All initial beams are ^/-polarized with transverse Gaus- 
sian distributions to guarantee the excitation of SPPs 
modes. The structure parameters are the same as in the 
MS simulation, with z G [0, 30] /im and perfectly matched 
layers at the boundaries. The analytical results without 
loss and the numerical results with loss are presented in 
Fig. 03 The data are shown for the x — z plane in the 
middle of the gap, y = ho/2 = 25 nm. 

For an on-axis input beam (xo = 0) [Fig. [3^a, b, e, f)], 
one can observe periodic beam focusing with a period of 

F = ir/n = 7ryjRoho(l + h b/a) . (6) 

The points of maximum intensity correspond to the fo- 
cusing when all excited modes interfere constructively. 
There is an initial phase shift depending on the ratio of 
the incident beam width and the characteristic width of 
the plasmonic waveguide wo/f]o- When these two widths 
match wo = 770, one can observe self-collimation effect - 
diffractionless propagation of the beam [see Fig.[3jc, g)]. 



For an off-axis input beam (xo 7^ 0) there are transverse 
oscillations [see Fig. [3jd, h)], and the plasmon polariton 
beam propagation exhibits an analog of optical pendulum 
behavior in a parabolic potential. The comparison with 
FDTD results suggests that losses do not affect the beam 
propagation except for attenuation along z direction. 

Eq. (|6]) has some important implications for propaga- 
tion of polychromatic SPPs beams. As a is practically 
wavelength independent while b varies slowly [Fig. [T^c)], 
the focusing strengths Q and the oscillation period F 
are nearly the same for different wavelengths. For ex- 
ample, in Fig. IHa) we show both the theoretical (solid 
line) and FDTD (circles) results for an input beam with 
wo — 800 nm > 770 for three different thicknesses ho in the 
wide wavelength range 0.5 — 1.0 /im. The different color 
beams focus at the planes of z = F\ = F/2. The dis- 
persionless characteristic of F enables broadband SPPs 
focusing and distinguishes our plasmonic structure from 
other highly wavelength dependent structures @, 0] . In 
Fig.[4](b-d) we show the field distribution |i£| on the x — z 
plane of y = ho/ 2 = 15 nm for the three points marked 
in Fig. Ufa) at ho = 30 nm. False colors are used to in- 
dicate schematically the three different wavelengths. It 
is important to note that the focal planes of all spec- 
tral components coincide within 500 nm (less than one 
wavelength deviation in the spectral range 0.5 — 1.0 /im) 
and are practically indistinguishable. Note that the focal 
depth is in the range of 2 — 3 /im for the different spectral 
components. Fig. [4] proves that this structure acts as a 
lens for broadband polychromatic light focusing, a func- 
tionality that has never before been discussed in other 



4 




-2 2 -2 2 
x (|iim) x (nm) 



FIG. 5: (Color online) (a,b) FDTD results for the field dis- 
tribution, \E\, in the x — z plane y = ho/2 (h = 40 nm and 
20 nm) and A = 0.7 /im. White dashed lines indicate the focal 
plane, (c) FWHM in the focal plane vs. ho. Marked points 
correspond to the cases (a) and (b), respectively. Insets: field 
distributions \E\ in the two marked focal planes. 



beam, bringing a potential functionality for manipulation 
of ultra-short optical signals. The focusing capability of 
this plasmonic lens could be further improved by tapering 
the dielectric layer along the propagation direction [20|. 
We anticipate the proposed structure is a promising can- 
didate for broadband plasmonic applications including 
subwavelength white light imaging, polychromatic plas- 
mon solitons, and ultrashort pulse plasmonic nanolasers. 
Furthermore, our ideas can be extended beyond the field 
of plasmonics, to include quantum particles in parabolic 
potentials, considering the similarity between Shrodinger 
equation and paraxial wave equation. 

We thank A. A. Sukhorukov, Z. Xu, A. S. Desyatnikov, 
C. G. Poulton, and A. E. Minovich for useful discussions, 
and acknowledge a support from the Australian Research 
Council and the NCI Merit Allocation Scheme. 



plasmonic structures. From Fig. [4] and Eq. (j6j), it is clear 
that one can control the position of the focusing points 
by varying the thicknesses ho and/or the wire curvature 
Ro, which enables wide reconfigurability of our structure. 

Finally, we investigate the focusing resolution of this 
structure, which is characterized by the full width at half 
maximum (FWHM) along x direction at the focal plane. 
(Along the vertical direction the resolution can be ap- 
proximately characterized by the thickness of the dielec- 
tric layer as shown in the insets of Fig. 0) We define the 
wavelength of the SPPs, A spp = A/n e ff, which is smaller 
than the light wavelength in the dielectric, X p = A/n^. 
What makes SPPs special is that its wavelength could 
be made as small as required, even vanishing for local- 
ized surface plasmons [2]. In our lens-like structure, if 
we want to improve the focusing resolution, we could 
decrease A spp . For a fixed vacuum wavelength, the sim- 
plest approach to increase n e s is to make the thickness ho 
smaller. Fig. [5] shows the dependency of FWHM along x 
at the focal plane for a fixed wavelength A = 0.7 jam and 
wo = 800 nm. Fig.JS^a, b) shows the field distribution \E\ 
along propagation for the two points marked in Fig. [5fc) 
in the x — z plane of y = ho/2 (20 nm and 10 nm, respec- 
tively). The white dashed lines indicate the focal planes. 
Insets (i) and (ii) show the transverse field distribution 
of \E\ at the two focusing planes. As shown in Fig. [5j 
the FWHM could be as small as 100 nm. This could 
be further improved by decreasing Ro and/or increasing 
the focusing strength Q. Our results clearly show that 
the light is fully confined in the dielectric along y and 
trapped in the optical parabolic potential along x. 

In conclusion, we have suggested the concept of poly- 
chromatic plasmonics and demonstrated a broadband 
plasmonic lens based on a parabolically modulated MDM 
structure. We have shown the plasmonic lens allows for 
complete three-dimensional subwavelength control of a 
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